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SUMMARY 


Radiance  spectra  for  two  hot -gas  samples  have  been  computed  using 
quasi -monochromatic  calculation  techniques  and  data  from  the  Air  Force 
Cambridge  Research  Laboratories  (AFCRL)  line  atlas  augmented  with  spectral 
data  on  HCl  and  HF.  Apparent  radiance  at  the  top  of  the  atmosphere  has  also 
been  computed  as  has  the  effective  average  atmospheric  transmittance,  which 
depends  on  the  source. 

The  adequacy  of  the  AFCRL  atlas  for  hot-gas  modeling  has  been  explored 
by  comparison  of  experimental  and  computed  spectra.  The  atlas  was  found  to 
contain  insufficient  hot  lines  for  accurate  modeling  of  the  wings  of  bands. 

Effort  should  be  directed  toward  providing  an  adequate  line  atlas  for  hot-gas 
modeling. 

The  HF  and  HCl  band  line  parameters  include  hot  lines,  but  omit  the 

2 2 
H isotopic  species.  Although  the  normal  concentration  of  H is  only  0.0156 

percent,  the  lines  are  strong  enough  to  ap  • •*  in  the  spectra.  These  lines 

should  be  added  in  the  future.  A more  app-  iate  line -width  parameter  for 

high  J values  is  also  needed. 

The  Voigt  line  shape  has  not  been  used  because  of  cost.  Efficient 
approximations  to  the  Voigt  function  should  be  obtained. 
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I.  INTRODUCTION 


High -re solution  spectra  of  the  radiation  from  a combustion-product 
source  viewed  through  an  atmospheric  absorption  path  are  required  for  several 
applications.  An  important  example  is  the  remote  sensing  of  missiles  or 
aircraft  from  their  exhaust  plume  radiation.  The  measurement  of  these 
radiance  spectra  would  require  the  use  of  research-laboratory -type  instru- 
ments  because  the  individual  spectral  lines  are  less  than  0.  1 cm”  wide. 
Because  measurements  are  difficult  in  the  field,  recourse  has  been  made  to 
computations.  If  the  pressure,  temperature,  and  species  concentrations  are 
sufficiently  well  defined  in  the  plume  and  atmosphere,  and  if  the  required 
spectral  coefficients  are  available,  the  radiance  can  be  computed  by  straight- 
forward but  tedious  means.  The  procedures  developed  to  make  these  compu- 
tations at  The  Aerospace  Corporation  were  described  in  an  earlier  report.  ^ 

The  spectral  coefficients  for  the  atmospheric  gases  were  obtained  from 

2 

the  Air  Force  Cambridge  Research  Laboratories  (AFCRL)  compilation. 
Although  this  line  compilation  was  never  intended  for  plume  modeling,  it  is 
the  most  comprehensive  compilation  available  and,  therefore,  was  used  in 
the  present  study.  The  compilation  has  been  revised  since  the  comparisons 
presented  in  Ref.  1.  In  Section  II,  the  revised  version  is  compared  with 
experimental  measurements. 

The  spectra  of  sources  that  contain  HCl  and  HF  were  investigated.  The 
spectral  coefficients  for  these  gases  were  calculated  from  basic  molecular 
parameters.  These  computations  and  a comparison  of  the  results  with  experi- 
mental line  parameters  and  experimental  spectra  are  given  in  Section  III. 

The  calculated  spectra  are  given  in  Section  IV. 


n.  AFCRL  LINE  ATLAS  COMPARISON 

The  AFCRL  compilation  of  spectroscopic  line  parameters  was  originally 
designed  to  study  the  propagation  of  infrared  radiation  through  the  atmosphere. 

The  atlas  is  adequate  for  these  studies  in  the  2.7-jim  regions,  as  the  compari- 
sons in  Ref.  1 indicated.  Since  the  publication  of  Ref.  1,  revised  editions  of 
the  line  atlas  have  been  received  from  AFCRL.  In  the  revised  version,  the 
error  in  the  010-130  water  band  strength  has  been  corrected,  several 
additional  water  bands  have  been  added  in  the  2.7-pim  regions,  and  the  CO^ 
parameters  have  been  significantly  reworked  throughout  the  infrared  spectral 
region.  Because  of  these  changes,  it  was  deemed  beneficial  to  repeat  the 
comparisons  made  in  the  previous  study  as  well  as  to  add  seme  additional 
comparisons. 

The  AFCRL  atlas,  as  revised,  appears  to  be  adequate  for  atmospheric 
paths.  The  comparisons  indicate  that,  in  spite  of  the  hot  bands  that  have 
been  added,  the  line  atlas  is  still  only  partially  adequate  for  modeling  gases 
with  temperatures  greater  than  1000  K. 

The  conditions  of  the  various  test  spectra  are  shown  in  Table  1.  Most 
of  the  comparisons  are  in  the  2. 7 -pm  spectral  region;  the  resiilts  of  the 
comparisons  are  discussed  below.  Tests  1 through  4 were  conducted  for  the 
same  conditions  as  the  tests  in  Ref.  1. 

A.  ATMOSPHERIC  CO^  PATH  (TEST  NO.  1)  { 

I I 

I Experimental  CO^  spectra  characteristic  of  atmospheric  paths  have 

i been  obtained  by  Burch,  ^ and  his  sample  10  was  chosen  for  comparison, 

j j The  solid  curve  in  Fig.  1 represents  the  experimental  results  and  was 

i ; obtained  by  differentiating  a tabulated  spectrum  from  Ref.  3;  the  dashed  curve 

i j indicates  the  result  obtained  in  the  earlier  study.  The  results  obtained  with 

• I the  most  recent  line  atlas  are  shown  as  the  broken  curve.  In  general,  the 

! revised  atlas  gives  results  that  approximate  those  of  the  old  atlas  and  are  in 

I reasonable  agreement  with  the  experimental  values.  Both  computed  spectra 

i I 

j 

^ I 

■ I 
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I III. 


Table  1.  Exper-mental  Conditions  for  Comparison  with  Computed  Spectra 


u 


FREQUENCY,  cm 


WAVELENGTH,  /xm 


Comparison  of  Computed  and  Experimental  Low- 
Resolution  CO2  Spectra.  The  solid  curve  is  the 
experimental  curve  of  Burch^  for  a homogeneous 
path  with  conditions  characteristic  of  an  atmo- 
spheric path.  The  dashed  curve  was  obtained 
from  the  AFCRL  line  atlas  used  by  Randall.  ^ 

The  1 roken  curve  was  computed  for  the  same 
conditions,  except  that  the  revised  AFCRL  line 
atlas  was  used. 


were  obtained  by  convolving  a 2.  5-cm  full -width -at -half -maximum  (FwHM) 
triangle  with  the  high-resolution  spectrum  computed  with  the  INHOM  program 
and  the  line  atlas  indicated. 

A highe r -resolution  spectrum  for  this  same  sample  is  .=’^own  in  Fig.  2. 
The  solid  experimental  curve  was  reproduced  from  Ref  3 by  procedures  out- 
lined by  Randall.  ^ The  computed  curves  were  o’  tained  by  convolving  a 
0.  8-cm"^-wide  FWHM  triangular  function  with  the  high-resolution  spectrum. 
The  dashed  curve  is  based  on  the  earlier  line  atlas.  The  results  obtained 
with  t'iie  revised  atlas  (broken  curve)  are  shown  only  in  the  two  regions  where 
they  differ  most  significantly  from  the  earlier  results.  Although  the  differences 
are  real,  they  are  not  large  enough  to  be  of  great  concern.  The  apparent 
eloser  agreement  of  the  early  atlas  results  with  the  experiment  results  is 
somewhat  surprising. 

B.  ATMOSPHERIC  H-O  PATH  (TEST  NO.  2) 


The  solid  curve  in  Fig.  3 represents  a portion  of  the  experimental  H-O 

4 

spectrum  obtained  by  Burch,  Gryvnak,  and  Patty  and  identified  as  sample  39. 
Computed  curves  were  obtained  by  convolving  a 0.  5-cm  -wide  FWHM  triangle 
with  the  computed  high- re  solution  spectrum.  The  strength  of  spectral  lines  of 
the  010-030  band  at  the  positions  indicated  by  the  vertical  bars  in  Fig.  3 were 
too  large  in  the  previous  version  of  the  line  atlas.  The  dashed  curve,  which 
was  obtained  with  the  revised  atlas,  indicates  that  this  problem  has  been 
corrected,  and  reasonably  good  agreement  exists  between  the  measured  and 
experimental  data. 

C.  HOT  CO^  PATHS  (TESTS  NO.  3 AND  5) 

The  solid  curve  in  Fig.  4 represents  the  experimental  radiance 
spectrum  for  Test  No.  3,  a hot,  high-pressure  sample  of  CO^.  The  spectrum 

and  experimental  conditions  were  obtained  from  the  work  of  Simmons,  Yamada, 

5 -1 

and  Arnold.  The  computational  results  were  obtained  by  convolving  a 5-cm 

wide  FWHM  triangular  instrument  function  with  the  high- re  solution  spectrum. 

The  dashed  curve  was  obtained  with  the  earlier  AFCRL  atlas.  The  broken 


FREQUENCY,  cm 


Figure  2.  Comparison  of  Computed  and  Experimental 
Medium  Resolution  CO2  Spectra.  The  same 
conditions  apply  as  for  Fig.  1,  except  at 
higher  resolution.  The  spectral  segments 
are  those  where  the  spectrum  computed 
from  the  revised  atlas  (broken  curve) 
deviates  most  significantly  from  the  pre- 
viously computed  (dashed  curve)  and 
experimental  results  (solid  curve) 
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EXPERIMENT 


FREQUENCY,  cm 


WAVELENGTH,  /^m 


Figure  4.  Comparison  of  Computed  and  Experimental  Hot  CO2  Spectrum. 

The  solid  curve  is  the  radiance  observed  by  Simmons^  from  a 
hot,  high-pressure  CO2  sample.  The  dashed  curve  was  com- 
puted for  the  same  conditions  but  with  the  line  atlas  used  in 
Ref.  1.  The  broken  curve  was  computed  from  the  revised 
line  atlas. 


curve  is  based  on  the  revised  atlas.  Tn  the  wings  of  the  band,  many  lines  are 
clearly  missing,  even  in  the  revised  atlas.  The  deviation  of  the  computed 
spectra  from  the  experimental  spectrum  in  the  center  of  the  band  indicates 
that  the  line  density  or  widths  or  both  are  not  great  enough  in  the  atlas  to 
make  the  sample  optically  thick,  whereas,  experimentally,  the  sample  is 
optically  thick. 

A hot  CO-  sample  (Test  No.  5)  for  which  the  path  is  not  optically 

^ 6 

thick  in  the  2.  7 -pm  band  was  studied  experimentally  by  Burch  and  Gryvnak. 

Their  result  is  shown  as  the  solid  curve  in  Fig.  5.  The  computational  result 

(dasl  jd  curve)  was  obtained  by  convolving  a 20 -cm  ^ -wide  FWHM  rectangular 

instrumental  function  with  a high -resolution  spectrum  obtained  with  the  INHOM 

program  for  the  conditions  listed  and  the  revised  line  atlas.  Again,  the  already 

noted  deficiencies  of  the  atlas  in  the  wings  of  the  CO^  band  and  the  inadequacies 

near  the  band  center  are  apparent. 

D.  HOT  H^O  PATH  (TEST  NO.  4) 

The  solid  curve  in  Fig,  6 represents  the  experimental  hot -water 
spectrum  obtained  from  the  work  of  Simmons,  Arnold,  and  Smith,  as 
described  by  Randall.  The  broken  crrve  is  the  computed  spectrum  obtained 
by  convolving  a 7.  5 -cm  ^ -wide  FWHM  triangular  instrumental  function  with 
the  high -re solution  spectrum  obtained  with  the  revised  line  atlas.  The  error 
caused  by  the  use  of  an  incorrect  030-010  band  intensity  in  the  previous  atlas 
has  been  corrected,  and  the  differences  between  the  computed  and  the  experi- 
mental spectra  are  now  in  the  expected  direction.  The  computed  transmittance 
is  higher  because  of  the  lack  of  the  hot  lines  in  the  atlas. 

Comparison  of  Figs.  5 and  6 indicates  that  the  AFCRL  atlas  is  better 
capable  of  predicting  the  transmittance  of  hot  H^O  than  hot  CO^.  Apparently, 
the  percentage  of  water  lines  with  significant  strength  at  high  temperatures 
included  in  the  atlas  is  greater  than  the  percentage  of  CO^  lines. 
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Figure  5 


3700 

FREQUENCY,  cm’^ 


Comparison  of  Computed  and  Observed  Hot  CO2 
Transmittcince.  The  solid  curve  is  the  spectrum 
observed  by  Burch.  ° The  dashed  curve  was 
computed  from  the  revised  line  atlas. 
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TRANSMITTANCE 


Figure  6.  Comparison  of  Computed  and  Observed  Hot  H2O  Transmittance. 

The  solid  curve  is  the  experimental  curve  of  Simmons  .T  The 
dashed  curve  was  computed  from  the  revised  line  atlas. 
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E.  HOT  CO^  PATH,  4.  3-|im  SPECTRAL  REGION 

(TEST  NO.  6) 

The  solid  experimental  curve  in  Fig.  7 was  obtained  at  the 
Environmental  Research  Institute  of  Michigan  (ERIM).  The  dashed  curve 
was  obtained  by  convolving  a 5-cm~^-wide  FWHM  triangular  instrument 
function  with  the  high  - resolution  spectrum  computed  with  the  revised  .^FCRL 
atlas  and  with  a pressure -broadened  line  shape.  For  the  pressures  involved 
in  this  experiment,  a Voigt  profile  is  more  appropriate;  however,  because 
the  Voigt  profile  is  more  complicated  mathematically,  and  hence  requires 
more  computation  time,  only  a portion  of  the  spectrum  was  computed  with 
the  Voigt  profile.  The  results  are  indicated  by  X's  in  Fig.  7.  Although  the 
use  of  a more  correct  profile  improves  the  agreement  with  experiment,  it  is 
clear  that  the  AFCRL  atlas  is  inadequate  for  accurate  plume  modeling  in  this 
spectral  region  (4.3-pm)  also. 

F.  NONISOTHERMAL  HOT  H^O  PATH  (TEST  NO.  7) 

The  spectra  for  the  only  nonisothermal  path  considered  in  evaluating 
the  AFCRL  line  atlas  are  shown  in  Fig.  8.  The  experimental  spectra  and  gas 
sample  conditions  were  obtained  from  the  work  of  Simmons,  Yamada,  and 

5 

Arnold.  The  temperature  profile  for  the  60 -cm -long  optical  path  is  shown 
in  Fig.  9.  The  open  data  points  denote  the  measured  temperatures  along  the 
path,  and  the  step -function  curve  is  the  profile  used  for  the  computations. 
Both  radiance  and  transmittance  spectra  are  shown  in  Fig.  8 inasmuch  as  the 
radiance  depends  on  temperature  through  the  blackbody  function  as  w'ell  as 
through  the  line  parameters.  The  theoretical  spectra  were  obtained  by  con- 
volving a 3.  5-cm~^ -wide  FWHM  triangular  instrument  function  with  the  high- 
resolution  spectrum  computed  with  the  revised  AFCRL  atlas  and  are  shown 
as  smooth  solid  curves  in  Fig.  8. 


F.  S.  Simmons,  private  communication  (9  August  197  1). 
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Comparison  of  Computed  and  Observed  Hot  CO2  b’pectra  in  the 
4.3-pm  Region.  The  solid  curve  was  obtained  experimentally 
at  ERIM.  The  dashed  curve  was  computed  from  the  revised 
AFCRL  atlas  with  use  of  a Lorentz  line  shape.  The  X -denoted 
curve  was  computed  from  the  same  line  atlas  with  use  of  the 
more  appropriate  Voigt  line  shape. 


FREQUENCY,  cm  ‘ 


Figure  8. 


Comparison  of  Experimental  and  Computed  Characteristics 
of  a Nonisothermal  H2O  Sample.  The  smooth  curves  are  the 
computed  transmittance  and  radiance  for  a path  with  the 
temperature  profile  shown  in  Fig.  9.  The  closed  data  points 
denote  comparable  experimental  points  from  the  work  of 
Simmons.  ^ 


DISTANCE,  cm 


Figure  9.  Temperature  Profile  for  Nonisothermzil  H2O  Computation. 

The  measured  temperatures  eilong  the  Scimple  cell  used  to 
obtain  the  experimental  points  in  Fig.  8 are  shown  (O). 
The  horizontal  line  segments  denote  the  temperature  pro 
file  that  was  used  to  simulate  the  sample  for  computation. 
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The  measured  spectra  in  Ref.  5 are  plotted  such  that  reproduction  j 

on  the  scale  of  Fig.  8 would  not  produce  a very  meaningful  comparison. 

Within  the  spectral  region  involved,  there  is  a one-to-one  c jrrespondence 
between  the  spectral  maxima  and  minima  in  the  measured  and  computed 
spectra.  These  extreme  values  were  obtained  from  Ref.  5 and  are  represented 
by  the  closed  data  points  in  Fig.  8.  The  solid  curves  connecting  these  points 
are  included  only  as  a visual  convenience  in  following  the  curve  and  do  net 
represent  intermediate  values  actually  measured. 

In  this  spectral  region,  the  agreement  between  experiment  and 
computation  is  generally  reasonable,  although  not  perfect.  In  the  trans- 
mittance plots,  the  tendency  is  for  the  computed  results  to  xinderevaluate 
the  absorption.  This  tendency  is  consistent  with  the  results  discussed  in 
Section  IID  above. 


III.  HCl  AND  HF  SPECTRAL  LINE  PARAMETERS 


m 


The  computation  of  spectral  r.  diance  from  plumes  that  contain  HF 
and  HCl  requires  that  the  line  position  v»  strength  S,  width  and  lower 
energy  level  E'^  be  available  for  all  lines  of  interest.  For  hot  plumes,  the 
lines  of  interest  include  not  only  the  transitions  that  involve  the  first  few 
rotational  levels  of  the  v = 1 to  v = 0 band,  which  have  been  studied  exten- 
sively in  the  laboratory,  but  also  transitions  that  involve  higher  vibrational 
and  rotational  levels.  Fortunately,  HF  and  HCl  are  simple  diatomic  mole- 
cules that  have  received  extensive  theoretical  analysis.  Procedures  are 
therefore  avilable  for  computing  the  energy  levels,  line  positions,  and 
strength  from  basic  molecular  constants.  Pressure  shifts  of  line  positions 
are  well  documented  for  these  gases,  but  were  ignored  in  the  present  study 
because  of  the  relatively  low  pressures  that  were  used.  The  theoretical 
prediction  of  the  pressure-broadening  parameter  is  more  difficult.  In  the 
present  study,  empirically  determined  line  widths  were  used,  where  avail- 
able, for  low-J  value  lines.  For  other  lines,  a constant  line  width  was  used. 

In  this  section,  the  basic  equations  used  to  determine  the  desired  line 
parameters  are  presented  together  with  the  values  of  the  molecular  constants 
employed.  Then,  the  accuracy  of  the  various  line  parameters  as  determined 
by  comparison  with  experimental  data  is  discussed,  and  a computed  HF 
spectrum  is  compared  with  an  experimental  spectrum. 

A.  BASIC  RELATIONS  AND  DATA 


1.  ENERGY  LEVELS 

The  energy  levels  or  term  values  for  a vibrating  and  rotating  diatomic 
molecule  have  been  studied  extensively.  The  theory  is  well  known  and  readily 

g 

available.  The  effects  of  the  interaction  of  vibration  and  rotation  are  accu- 
rately included  in  the  Dunham  formulation,  in  which  the  term  value  0 _ is 


HSCEDIbXJ  PrtOS  BLANK-HOT  ?IL‘tSD 
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given  in  terms  of  a double  series  in  the  vibration  quantum  number  v and 
the  rotational  quantum  number  J 


i=L  m=M 

®V,J  Z (1) 

i=0  m=0 


The  known  as  Dunham  coefficients,  are  usually  determined  by  data- 

fitting  techniques  from  the  positions  of  spectral  lines  in  carefully  calibrated 
spectra.  The  values  used  in  the  piesent  study  for  HF  and  HCl  are  given  in 

O c 

Table  2.  Sets  of  coefficients  for  HF  and  HCl  are  available  from  the 

‘ 9 i 0 ^7 

literature.  ’ A complete  experimental  analysis  for  HCl^  ' was  not  found. 

Therefore,  the  HCl  set  of  coefficients  was  derived  fron-  the  HCF^  by 

means  of  the  Dunham  isotopic  equation 


y]  = Y, 
i,m  i,m  ' 


where  p is  the  reduced  mass,  and  superscript  i implies  isotope  and  is  not 
an  exponent.  The  appropriate  values  of  q = i + 2m  are  given  in  Table  2. 
The  are  of  order  unity.  For  HCl^^  and  HCl^^ 


(B  /'D  (p  - P^)/p^  « - 2 X 10 

C 6 


Consequently,  the  terms  involving  this  ratio  were  ignored.  The  value 
(p/p  ) = 0.99924  was  used  for  deriving  HCl  coefficients  from  HCl  ^ 

values . 

These  term  values  include  the  zero-point  energy  of  the  oscillator. 

The  approximation  for  the  partition  function  in  the  present  programs  refers 
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Table  2.  Dunham  Coefficients  for  HF  and  HCl 


Coefficient 

HF,  ^ cm  ^ 

„„,35  b -1 

HCl  , cm 

q 

37 

HCl 

Value,  cm  ^ 

4138. 73 

2990.9463 

1 

2988. 6732 

^20 

-90.05 

-52.81856 

2 

-52.73831 

^^30 

0.  932 

0.22437 

3 

0.22386 

^40 

-1.42  X 10"^ 

-0.01218 

4 

-0.01214 

"^50 

-5.9  X lo”"^ 

5 

^01 

20.9555 

10.593416 

2 

10. 577320 

^11 

-0.7958 

-0.307181 

3 

-0. 306481 

^21 

1.  182  X 10“^ 

1.7724  X 10"^ 

4 

1.7670  X 10"^ 

^31 

-3.  11  X lO"'^ 

-1.201  X lO"'^ 

5 

-1.  196  X 10"'^ 

^41 

-5.8  X 10“^ 

6 

"^02 

-2.  153  X 10'^ 

-5.31936  X lO"'^ 

4 

-5.  30321  X lO"'^ 

^12 

6. 23  X 10"^ 

7.510  X 10"^ 

5 

7.482  X 10"^ 

^22 

-2.06  X 10'^ 

-4.00  X 10“^ 

6 

-3.98  X 10"^ 

^03 

1.68  X 10"^ 

1.74  X 10"® 

6 

1. 732  X 10"^ 

-6.  5 X 10“^ 

6.34  X 10“^° 

7 

6.31  X 10"^° 

^04 

-1.0  X 10'^^ 

-9.93  X 10"^^ 

8 

-9. 87  X 10"^^ 

Ref.  9 
^Ref.  10 


all  energy  levels  to  the  lowest  vibrational  level.  Therefore,  the  lower 
energy  level  required  from  the  line  atlas  for  population  computation  is 


E''  = 0 T - n 

v,J  v,J  0,0 


2.  LINE  POSITIONS 

The  frequency  of  a spectral  line  is  given  by  the  difference  in  term 
values  for  the  upper  and  lower  energy  levels  between  which  the  transition 
occurs.  The  vibrational  quantum  number  v can  change  by  any  integral  value. 
Vibrational  bands  with  Av  > 2 lie  outside  the  spectral  region  of  present 
interest  and  are  excluded  from  the  present  study,  except  for  the  v = 0 to 
V = 2 band.  Because  of  anharmonicities  in  the  moleciilar  vibration,  the 
Av  = 1 bands  are  not  all  superimposed  on  each  other.  Four  Av  = 1 bands 
with  lower-level  v numbers  equal  to  0,  1,  2,  3 are  included  in  the  present 
study.  For  the  molecules  of  interest  here,  the  rotational  quantum  number  J 
can  change  by  ±1.  The  R branch  of  the  vibrational  band  corresponds  to 
AJ  = +1  when  a photon  is  absorbed:  the  P branch  corresponds  to  AJ  = -1 
when  a photon  is  absorbed.  The  conventional  notation  is  to  indicate  AJ  by  a 
letter  and  the  J value  of  the  lower  state.  With  this  notation,  the  frequencies 
of  the  lines  may  then  be  specified 


v",J  ~ v',J+l  ” v",J 


'^v',v",J  ""  ®vLJ-l  ■ ®v'LJ 


(R  Branch) 


(P  Branch) 


Note  that  lore  is  no  P(0)  line. 

LINE  STRENGTHS 

Line  strengths  are  defined  in  several  ways.  For  clarity,  we  first 
carefully  define  how  line  strengths  are  used  in  the  computations  performed 


by  the  program  described  in  Ref.  1.  We  then  discuss  how  these  line 
strengths  are  computed  from  basic  molecular  constants. 

A homogeneous  gas  sample  at  temperature  T with  a single  absorption 
line  at  Vq  will  have,  for  radiation  at  frequency  v,  a transmittance 

T = exp  t-S(T)  Wf  (v,Vq)]  (4 

The  line  shape  function  is  denoted  by  f (v»  Vq)  and  will  have  the  unit  of 
reciprocal  frequency.  For  a pressure-broadened  line  of  width  a. 

f (v,Vq)  = a/  jir  j^a^  + (v  - 

W denotes  a measure  of  the  number  of  absorbing  molecules  present  and  is 
a column  density  with  unit  of  molecules  per  unit  area.  The  desired  unit  for 
the  line  strength  is,  then,  frequency  per  unit  column  density.  Specifically, 
for  the  present  calculation,  S has  the  unit  cm  /{molecules  cm  ). 

In  the  study  of  gases,  where  the  ideal  gas  law  is  valid,  the  column 
density  W can  be  related  to  the  temperature  T,  pressure  P,  and  length  L 
of  the  sample  by  means  of  the  Boltzmann  constant  k 

W = PL/(kT)  (6 

\ 

If  this  is  combined  with  Eq.  (4) 

T = exp  [-S(T)  PL  f(v,VQ)/(kT)]  (7 

The  ratio  <y  = S(T)/kT  is  also  called  the  line  strength  by  many  a\’.thors;  the 
usual  unit  is  atm  ^ cm  ^ 


The  relationship  between  these  two  definitions  of  line  strength  is 


a = C S/T 
o 

C = 7.339083  X 10^^  molecules  K cm  ^ atm 
o 


Except  for  comparison  with  experimental  results,  the  line  strength  S is 
exclusively  in  this  report. 


The  line  strength  S depends  on  the  fraction  of  the  molecules  in  the 

energy  level  from  which  the  transition  begins  and  the  matrix  element 
2 

R »T/  urn  by  the  relations 
v'J  v"J''  ^ 


v'v"J'J"  " 3hc 

^v'J'v"j"  ''O 

(-J 

for  P branch  transitions 

m = 

fj+1 

for  R branch  transitions 

""  ''v',v",  J',  J" 
g(vQ  T)  = 1 - exp  (-  hcvg/kT) 


Z =2^  (2j+l)exp(-9^  hc/kT) 

v,J 


^ i 
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(8) 

used 

(9) 

(10) 

(11) 

(12) 

(13) 


where  h is  Planck's  constant,  6.626196  X 10  erg  sec;  c is  the  velocity 
of  light,  2.997925  X 10^^  cm  sec  S and  k is  Boltzmann's  constant, 

1.380622  X 10"^^  erg  K"  ^ 

The  line  strength  S,  as  defined  by  Eq.  (9),  is  a complicated  function 
of  temperature.  Explicit  account  of  all  the  various  temperature  variations 
would  be  time  consuming.  Therefore,  for  the  spectral  computations,  the 
line  strength  S®  computed  at  the  standard  temperature  Tq  = 296  K by  means 
of  Eq.  (9)  is  scaled  to  other  temperatures  by  the  relation 


S = S"  g(v,T)  exp 


hcE" 

k 


(1/Tq  - 1/T)  /z(T) 


where  z(T)  is  the  ratio  of  the  partition  function  at  T to  the  value  of  the 
partition  function  at  T^.  In  order  to  avoid  the  summation  implied  by  Eq.  (13), 
z is  approximated  by  the  product  of  separate  ratios  for  the  v-Tirational  and 
rotational  partition  functions 

z(T)  = Q„  (T)  Q (T)  = (T/TJ  [1  - exp  (-  Y,„hc/kT)]  (15) 


where  Y,_  is  used  as  the  vibrational  frequency.  No  difference  between 
35  37 

HCl  and  HCl  is  recognized  in  calculating  z(T). 

Except  for  the  matrix  elements,  all  the  factors  in  Eq.  (9)  can  be 
computed  from  the  term  values  defined  previously.  Computation  of  the 
matrix  elements  requires  computation  of  the  integral  of  the  radial  wave 
functions  ^(r)  for  the  initial  and  final  states  with  the  electric  dipole  moment 
function  M(r) 
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M(r)  is  usually  expanded  in  a power  series  about  the  equilibrium  nuclear 
separation  r 


M(r)  M.  (r  - r 


Generally,  the  0(r)  functions  can  be  generated  m.ost  accurately  by  numerical 

solutions  of  the  radial  Schroed’nger  equation.  The  integral  to  obtain  the 

dipole  matrix  elements  is  then  also  performed  numerically.  This  procedure 

35 

has  been  applied  to  several  molecules,  including  HF  and  HCl  , to  obtain 

12-  14 

accurate  line  strengths.  This  approach  requires  a significant  com- 

putational effort  independent  of  the  spectral  computation  of  primary  interest 

in  the  present  study.  Therefore,  we  elected  to  use  the  analytical  approxi- 

15 

mations  to  Eq.  (16)  formulated  by  Herman  and  Wallis. 

In  these  analytical  approximations  for  the  matrix  elements,  the 
maxtrix  element  is  written  as  a product  of  an  anharmonic  oscillator  factor 
D and  a coxrection  factor  F,  which  accounts  for  vibration-rotation 
interaction. 


v',v",j',J 


^ - "n  ^ F* 

" “ ■^v'^v'  v''v' ) 
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Formulas  for  the  D . are  available  in  terms  of  the  coefficient  of  the 

v^v'' 

first  power  in  the  dipole  moment  function  Mj  and  spectroscopic  constants 
defined  in  Eqs.  (19) 


D(3i  = Mj/(2a)' 


Di2  = Mj/cv 


0^3  = Mj/(2ff/3)^^^ 


D34  = Mj/(<y/2) 


°02  = Mjb/(2<y)' 


where 


Of  = u)  / (2B  r'^) 
e e e 


b = -0.5V"'  {1  + a^/(3YBg) 


Y = 2B  /(I) 
e e 


The  equilibrium  spacing  of  the  atoms  is  r ; the  equilibrium  rotation  constant 


B the  equilibrium  vibration  frequency  lu  Y.^;  and  the  rate  of 

G V A i J ^ ^ 

change  of  B with  v is  a ^ - Y, . . Note  that  a and  » are  different 
e e n e 


quantities. 


Formulas  have  also  been  computed  for  the  vib.'ation- rotation  inter- 
action correction  terms 


{Fq  = 1 - 26Ym  [^1  + 7.5  bY^^^  + 3Y  - 2.5  bY^^^/9 

+ 0.75  bY^^^  (m-  l)/0  + 0.375Y  (m  - l)/ej 
(Ff  = 1 - ^«Ym  jl  + 15  bY^^^  + 6Y  - 23  bY^^“/(40) 


+ 0.75  bY  ' (m-  l)/0  + 0.375Y  (m 


(cont. ) 
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(F^  3)^^^  = 1 - 20Ym  - 450bY^^^  m - 180Y^m  - 0.75Y^m(m-  3) 
- i.5m^  bY^^^  + 39  mbY^^^/2 
(F^  4)^^^  = 1 - 20Ym  - 240Y^m  - 299embY^'^^/4 

+ 6lmbY^^^/2  - 0.75m^bY^^^  - 0.  75Y^  m (m  - 2) 
(Fq  1 (^/b)  [(1  - 0.750  + O.5b0/Y^^^)  - 0Ym 

+ 0.5  (2Y^^^/b  + 15/4)Y*^^  b(m  + 5) 


+ (3/16)Y^^^  b(2m-  3)j 


All  the  factors  except  0 = M^/M^r^  have  been  previously  defined. 

The  values  used  with  Eqs.  (9)  through  (13)  and  (18)  through  (21)  to 

obtain  the  line  strengths  for  the  present  work  are  given  in  Table  3.  The  data 

were  taken  from  the  references  given  in  the  table.  The  values  of  the  dipole 

37  35 

function  expansion  coefficients  for  HCl  were  obtained  from  the  HCl 

37 

values  by  multiplying  by  the  square  root  of  the  1:3  abundance  ratio  for  Cl 
35 

compared  with  Cl 
4.  LINE  WIDTHS 

Pressure  broadening  of  spectral  lines  results  from  collisions  between 
molecules.  The  relative  effectiveness  of  a collision  in  broadening  a line 
depends  on  the  molecules  involved  in  the  collision  and  on  the  energy  state  of 
the  molecules.  Collisions  between  molecules  of  the  same  type  and  in 
similar  rotational  states  are  the  most  effective  collisions  for  the  line 
broadening.  The  transition  involved  is  much  less  important.  These 
featores  are  illustrated  in  Figs.  10  and  11. 

The  general  similarity  of  the  dependence  of  line  width  on  the  rotational 
quantum  number  factor  |m|  (Eq.  (10)]  for  a large  variety  of  HCl  transitions 
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LINE-WIDTH  PARAMETER 


Figure  10.  Line-Width  Parameter  for  HCl  as  a Function  of  |m|.  The  line 
connects  values  used  for  the  line  parameters  obtained  in  the 
present  study.  Excep  for  □,  the  symbols  represent  experi- 
mental self-broadening  values  obtained  for  different  bands  and 
isotopes  by  various  investigators.  O,  HCl^^,  0-1,  R branch 
(Ref.  18);  •,  HCl^^,  0-1,  R branch  (Ref.  18);  A,  0-1  P branch 
HCi35  (Ref.  18);  A,  0-1,  P branch  HCP’^  (Ref.  18);  CD,  0-1 
P branch,  both  isotopes  (Ref.  21),  X,  0-2  band,  both  branches 
(Ref.  20);  +,  0-4  band  (Ref.  19);  □,  N2  broadening  of  HCl 
(Ref.  18). 


PARAMETER, 


o 

200  K 

o 

300  K 

□ 

513  K 

A 

1283  K 

Figure  11.  Measured  Line  Widths  of  the  R Branch  of  the  HCl  Overtone 
for  Several  Temperatures.  The  solid  curves  were  calcu- 
lated from  the  Anderson  theory  for  resonant  dipoles  with 
temperature -dependent  hard  cores.  Taken  from  Goldring 
and  Benesch. 





is  illustrated  in  Fig.  10.  The  open  and  closed  triangles  and  circles  from 

1 3 

the  work  of  Benedict  et  al.  represent  the  P and  R branches  of  the  1-0  band 
35  37 

of  HCl  and  HCl  ; the  crosses  represent  the  4-0  band  of  HCl  as  reported 

19 

by  Benedict  et  al.  from  the  work  of  Lindholm;  the  X's  represent  the  2-0 

band  of  HCl  as  measured  by  Goldring  and  Benesch;^^  and  the  ellipses 

represent  the  P branch  of  the  1-0  band  as  measured  by  Babrov,  Ameer, 

2 1 

and  Benesch.  These  are  all  self-broadening  coefficients,  i.  e.,  coeffi- 
cients for  collisions  of  HCl  with  HCl.  The  solid  curve  joins  the  values 

selected  for  HCl  in  the  present  study.  N,  is  a much  less  efficient  broaden- 

^ 18 

ing  agent,  as  is  noted  by  the  data  taken  from  the  work  of  Benedict  et  al. 
for  nitrogen  broadening  and  denoted  by  open  squares  in  Fig.  10. 

The  shape  of  the  relationship  between  jml  and  broadening  is  tempera- 
ture dependent,  as  is  shown  in  Fig.  11,  which  was  taken  from  Ref.  20. 

The  broadening  parameters  for  the  R branch  of  the  HCl  fundamental  are 
plotted  in  the  figure  as  a function  of  jmj  for  a variety  of  temperatures.  The 
flattening  of  the  curve  as  the  temperature  is  increased  may  be  easily  under- 
stood in  a qualitative  sense.  The  most  effective  broadening  collisions  are 
between  molecules  in  similar  J states.  At  low  temperatures,  most  of  the 
molecules  are  in  a low  J state.  Thus,  the  probability  oi  collisions  with  a 
molecule  in  a similar  J state  is  high.  As  the  temperature  increases,  the 
distribution  in  J states  becomes  less  peaked,  and  any  given  collision  is  less 
likely  to  be  with  a molecule  in  a similar  J state. 

The  qualitative  features  illustrated  in  Fig.  10  guided  the  empirical 
assignment  of  line  width  parameters  in  the  present  study.  The  variation  in 
the  relative  shape  of  the  broadening  parameter  curve  with  temperature  was 
not  included  because  of  complications  in  using  this  information  in  actual 
spectral  calculations. 

The  line  width  parameters  assigned  for  HCl  and  HF  taken  from 
1 8 ZZ 

Benedict  et  al.  and  Smith  are  given  in  Table  4.  Line  widths  for  Iml 
values  greater  than  those  for  which  data  are  available  were  assigned  the 
constant  values  shown  in  the  table. 
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Table  4.  Line-Width  Parameters  for  HF  and  HCl 


HCl/ 

hf/ 

-1  -1 

-1,-1 

!ml 

cm  atm 

cm  atm 

.218 
.223 
.249 
.246 
.227 
. 186 
. 153 
. 127 
. 106 
0.0896 
0.0842 
0.0834 
0.066 
0.C6 


0.56 
0.67 
0.63 
0.50 
0.33 
0.24 
0. 15 
0. 10 
0. 10 


The  spectral  computation  program  computes  a composite  line  width 
a appropriate  for  a particular  path  segment.  The  computation  is  based  on 
the  partial  pressures  P.  of  the  various  gases  present  in  that  segment  and 
the  line  width  parameter 

a = ao  (VT)‘'2J^T1.P. 

i 

The  coefficients  T).  are  the  ratios  of  the  foreign- gas  broadening  to  self- 
broadening coefficients  for  the  various  gases  included  in  the  program.  The 

21 

values  of  T].  were  estimated  from  the  work  of  Babrov,  Ameer,  and  Benesch 

for  HCl  (Table  5).  The  values  of  T)  used  for  CO2  and  CO  are  the  average 

ratios  shown  in  Table  5.  All  gases  other  than  the  halide,  CO,  and  CO^  are 

assumed  to  be  air.  For  air,  T)  = 0.279  is  obtained  by  taking  a weighted 

average  of  the  valu.rs  for  N_  and  O..  in  Table  5.  Based  on  the  empirical  results 

23  ^ ^ 

of  Shaw  and  Lovell,  the  Tj  coefficients  derived  for  HCl  were  also  used  for  HF. 

B.  COMPUTED  LINE  PARAMETER  ACCURACY 

i"*  I 

,,i  The  calculated  HF  and  HCl  spectra  in  this  report  are  based  on  a set 

of  spectral  line  parameters  that  were  computed  from  data  presented  in  the 

preceding  section  for  all  J values  up  to  50  for  the  0-1,  1-2,  2-3,  3-4,  and 

0-2  bands  of  HCl  , HCl  ' , and  HF.  This  set  includes  all  lines  of  these 

bands  for  which  a 0.  5-m  path  through  a 1-atm  3000 -K  sample  of  the  gas 

would  have  a peak  absorptance  of  0.01  or  more.  With  an  arbitrary  cutoff 

in  J,  a number  of  weaker  lines  are  also  included.  A more  rational  selection 

criterion  was  applied  in  selecting  the  line  parameters  given  in  the  Appendix. 

There,  the  parameters  are  given  only  for  lines  that  yield  an  absorptance 
-2 

greater  than  10  for  a 1-m  path  through  a 3000-K  sample  of  the  gas  at  1-atm 
pressure. 
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Table  5.  Self-Broadening  Parameter  Derivation 


Line 

HCl -HCl  Line- 
Width  Parameter 
cm“^  atm"^ 

Ratio  Foreign-gas  Broadening/ 
HCl  Broadening 

CO2 

CO 

^2 

°2 

PI 

0.  207 

0.807 

0.565 

0.447 

0.341 

P2 

0.  233 

0.600 

0.429 

0.365 

0.218 

P3 

0.231 

0.445 

0.415 

0.313 

0. 166 

P4 

0.  227 

0.392 

0.368 

0.283 

0.  147 

P5 

0.  207 

0.351 

0.  327 

0.286 

0.  142 

P6 

0.  171 

0.385 

0.348 

0.254 

0.158 

P7 

0.  145 

0.346 

0.394 

0.250 

0. 150 

P8 

0.  124 

0.380 

0.280 

0.215 

0,173 

Average  value 

0.463 

0.391 

0.302 

0. 187 

The  accuracy  of  line  positions  is  about  0.04  cm  ^ for  HF  and  about 
0.005  cm”  ^ for  the  H7-1  positions.  The  lower  energy  levels  are  determined 
with  similar  precision.  The  line  strengths  are  accurate  to  within  20  per- 
cent for  low  V and  J values.  The  uncertainty  increases  for  higher  J and  v 
values.  The  pressure-broadening  coefficients  are  within  50  percent  for  low 
temperature  and  m < 10;  they  are  less  accurate  for  other  conditions.  These 
statements  of  accuracy  are  demonstrated  in  the  next  several  paragraphs  by 
comparison  of  calculated  line  parameters  with  experimental  values. 

The  line  frequencies  of  HCl  are  the  most  accurate  quantities.  The 

35 

Dunham  parameters  from  Ref.  10  were  obtained  for  HCl  . Our  calculated 

35  - 1 

HCl  line  positions  agree,  to  within  0.001  cm  , with  the  calculated  values 
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of  Rank  et  al.,  which,  in  turn,  agree  with  the  observed  frequencies  to  within 

- 1 35 

about  0.004  cm  . An  independent  check  on  both  the  HCl  Dunham  coeffi- 

37 

cients  and  the  conversion  to  HCl  coefficients  is  provided  by  the  data  in 

37  24 

Table  6,  where  observed  1-0  band  HCl  line  frequencies  are  compared 

with  the  values  calculated  with  the  parameters  given  in  Table  2.  The 

35 

differences  are  somewhat  larger  than  with  the  HCl  observations,  as  would 
be  expected,  but  the  accuracy  in  this  band  is  still  better  than  0.01  cm 

9 

The  accuracy  of  the  HF  observations  of  Mann  et  al.  is  lower  than 
that  of  the  HCl  observations.  Consequently,  the  Dunahm  coefficients  are 
given  with  lower  accuracy.  Comparison  of  computed  and  calculated  values 
indicates  an  accuracy  of  about  0.04  cm 

The  lower  energy  levels  are  based  on  the  same  data  as  the  line 
frequencies  and  hence  should  have  a comparable  accuracy.  This  is  several 
orders  of  magnitude  more  accuracy  than  is  required  for  scaling  the  very 
much  less  precise  line  strengths  with  temperabire. 

There  are  only  a limited  number  of  experimental  line  strength  values 

available  for  HCl,  and  these  are  restricted  to  the  lower  rotational  states. 

In  Fig.  12,  the  calculated  values  for  the  0-1  band  (solid  curve)  are  compared 

2 1 

with  the  experimental  values  of  Babrov,  Ameer,  and  Benesch  (squares) 

1 8 

and  the  experimenir!  values  of  Benedict  et  al.  (circles).  In  Fig.  13,  the 

calculated  values  for  the  0-2  band  (solid  curve)  are  compared  with  the 

25 

experimental  values  of  Benedict  et  al.  These  two  sets  of  experimental 
data  differ  from  each  other  and  from  the  calculated  values  by  10  to  20  per- 
cent. This  is  the  uncertainty  we  associate  with  the  line  strength  for  these 
bands.  For  higher  vibrational  states,  the  uncertainty  in  band  strength  can 
be  expected  to  be  somewhat  larger. 

Line  strength  at  high  temperature  can  be  expected  to  have  a decreased 
accuracy  because  of  the  scaling  approximations,  Eqs.  (14)  and  (15),  as  well 
as  the  extrapolation  of  the  basic  theory  to  higher  excitation  levels.  In 
Fig.  14,  line  strengths  for  the  v = 2 to  v = 3 band  of  HCl  are  compared  with 
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Figure  12.  Comparison  of  Computed  and  Experi- 
mental HCl  Line  Strengths  for  the 
0-1  Band.  The  values  obtained  in  the 
present  study  are  connected  by  a 
solid  curve.  Open  symbols  denote 
experimental  HCl^^  values;  closed 
symbols  represent  HCl^'  values. 
Values  denoted  by  circles  are  from 
Ref.  18,  those  denoted  by  squares 
are  from  Ref.  21. 
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Figure  13.  Comparison  of  Computed  and  Experi- 
mental HCl  Line  Strengths  for  the  0-2 
Band.  The  values  obtained  in  the  pre- 
sent study  are  connected  by  a solid 
curve.  Open  symbols  denote  HCl^^ 
experimental  values;  closed  symbols 
denote  HCl37  values.  All  experimental 
values  are  from  Ref.  25. 


Figure  14.  Comparison  of  Computed  and  Experimental  HCl 
Line  Strengths  for  the  2-3  Band.  The  values 
obtained  in  the  present  study,  scaled  to  1200  K 
as  discussed  in  the  text  and  reduced  to  units  of 
a at  300  K,  are  connected  by  a solid  curve.  Th 
experimental  values  in  these  same  units  re- 
ported in  Ref.  25  for  HCl^^(O)  and  KCl^'^  (•) 
are  also  shown. 


25 

the  values  obtained  by  Benedict  et  al.  at  1200  K and  values  scaled  from 
the  296  K line  strengths  by  use  of  Eqs.  (14)  and  (15).  Although  the  agree- 
ment is  poorer  than  for  the  lower  temperatures  and  vibrational  quantum 
nunribers,  all  of  the  experimental  points  are  still  within  50  percent  of  the 
calctilated  values. 

In  Fig.  15,  the  computed  line  strengths  for  the  0-1  band  of  HF  are 

Z6 

compared  with  experimental  measurements  of  Lovell  and  Herget,  Shaw 
23  27 

and  Lovell,  and  Kuipers.  The  comparisons  are  made  for  a temperature 
of  390  K because  most  measurements  are  made  in  this  temperature  range  in 
order  to  avoid  dimerization  of  the  HF.  The  computed  line  strengths  for 
296  K have  been  scaled  to  390  K by  means  of  Eqs.  (14)  and  (15).  Except  for 
the  P(8)  and  P(9)  lines,  the  strengths  all  agree  with  each  other  and  with  the 
computed  values  to  within  20  percent. 

In  Fig.  16,  the  computed  line  strengths  for  the  0-2  band  of  HF  are 

28 

compared  with  the  observations  of  Meredith.  hi  view  of  the  good  agree- 
ment in  Fig.  15,  the  rather  consistent  deviation  of  the  measured  and  com- 
puted values  for  the  low  jm]  values  is  puzzling,  and  we  are  unable  to  explain 
it.  However,  the  measured  and  calculated  values  are  still  within  30  percent 
of  each  other. 

Data  OR  iine  widths  for  HCl  have  been  summarized  in  Fig.  10.  The 
line -broadening  parameters  for  low  tempo xatures  and  values  of  m < 13  are 
probably  within  20  percent  of  the  true  values.  For  higher  rotational  values 
and  higher  temperatures,  greater  deviations  are  expected  between  the 
parameters  used  for  calculations  of  spectra  and  the  true  values. 

HF  line-width  parameters  are  usually  measured  at  390  K or  greater 

in  order  to  avoid  problems  with  HF  dimers.  The  extreme  values  observed 

22  23 

at  this  temperature  in  five  different  studies  by  Smith,  Shaw  and  Lovell, 

26  2T  28 

Lovell  and  Herget,  Kuiper,  and  Meredith  are  indicated  in  Fig.  17 

by  vertical  bars  and  listed  in  Table  7.  The  values  obtained  from  the  302  K 
22 

data  of  Smith  and  used  in  the  line-parameter  compilations  have  been  scaled 
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Figure  15.  Comparison  of  Computed  and  Experimental 
HF  Line  Strengths  for  the  0-1  Band.  The 
values  obtained  in  the  present  study  are 
connected  by  a solid  curve.  Experimental 
values  are  denoted  by  □,  Kuipers;^^  A,  _ 
Shaw  and  Lovell O,  Lovell  and  Herget. 
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Figure  16.  Comparison  cf  Computed  and  Experimental 
HF  Line  Strengths  for  the  0-2  Band,  The 
values  obtained  in  the  present  study  are 
connected  by  a solid  curve.  Experimental 
values  from  the  work  of  Meredith^®  are 
shown  (O). 
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Figure  17,  Summary  of  HF  Experimental  Line -Width  Parameter  Data. 

The  ver^lcal  bars  indicate  the  range  of  the  values  observed 
by  various  investigators  and  listed  in  Table  7.  The  values 
used  in  the  present  study  are  connected  by  straight  line 
segments, 
ii 


HF  Line-Width  Parameter  Data 


to  390  K as  in  the  program  of  Randall  and  are  shown  connected  by  straight 
line  segments.  As  with  HCl,  a 20-percent  variation  in  experimental  values 
is  observed.  The  effect  of  the  change  with  temperature  of  the  |ml 
dependence  of  the  line-width  parameter,  which  is  not  included  in  the  present 
calculations,  is  also  apparent. 

Of  the  several  parameters  in  the  halide  line-parameter  compilations, 

the  line  width  is  the  least  accurate.  Theoretical  models  of  pressure 

29 

broadening  are  available.  Although  the  results  are  not  in  complete  agree- 
ment with  measured  values  {Fig.  11),  the  application  of  a consistent  theory 
for  extrapolation,  rather  than  the  ad  hoc  procedures  used  here,  is  certainly 
preferable  for  improvement  of  the  line  width  values.  Simple  methods  for 
incorporating  the  variation  in  J dependence  with  temperature  should  also 
be  incorporated  in  the  spectral  calculation  programs. 

The  results  of  an  independent  study  similar  to  the  present  one  for 

30 

the  Av  = 1 bands  of  HF  and  DF  were  recently  published  by  Goldman  et  al. 
Their  basic  data  sources  were  the  same  as  those  used  here.  Their  line 
parameters  were  published  as  a table  of  line  strengths  (o)  at  1273  K.  The 
line  strengths  for  the  2-3  band  of  HF  are  compared  in  Fig.  18,  where  the 
solid  curve  represents  the  results  of  the  present  study,  and  the  dashed 
curve  represents  those  of  Goldman.  The  P-branch  results  are  indistin- 
guishable on  the  scale  of  the  figure.  Goldman’s  R-branch  results  are  a 
few  percent  lower  than  our  results,  but  certainly  agree  to  within  the  accu- 
racy quoted  for  either  set  of  results. 

C.  HF  SPECTRAL  COMPARISON 

The  final  measure  of  the  adequacy  of  a line-parameter  compilation 
for  spectral  computation  is  the  comparison  of  computed  and  experimental 
spectra.  In  Fig.  19  are  shown  the  spectra  computed  for  a 60-cm  sample 
of  HF  at  1 atm  and  with  the  inhomogeneous  temperature  distribution  shown 
in  Fig.  20.  The  spectra  presented  are  the  result  of  convolving  a 4.5-cm~  - 
wide  FWHM  triangular  instrument  function  with  the  computed  high- re  solution 


Figure  18,  Comparison  of  the  Computed  HF  Line  Strengths  of  Goldman, 
et  al.  at  1273  K with  Values  from  the  Present  Study  in  the 
2-3  Band.  The  solid  curve  represents  results  of  the  pre- 
sent study.  Goldman's  results  are  denoted  by  a broken 
curve,  which  is  indistinguishable  from  the  solid  curve  in 
the  P branch. 
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Figure  20.  Temperature  Profile  for  Nonisothermal  HF  Computation. 

The  measured  temperatures  along  the  sample  cell  used  to 
obtain  the  experimental  points  in  Fig.  19  are  shown  (O). 
The  horizontal  line  segments  denote  the  temperature  pro- 
file used  to  simulate  the  sample  for  computation. 


I 

spectrum.  These  conditions  agree  with  experimental  spectra  obtained  by- 
Simmons,  Yamada,  and  Arnold.^  Because  the  experimental  spectra  do 
not  lend  themselves  to  reproduction,  we  have  extracted  the  extreme  values 
from  the  experimental  spectra  and  plotted  them  in  Fig.  19.  The  agreement 
is  quite  satisfactory. 


IV.  HOT-GAS  SPECTRA 


In  this  section,  the  spectral  line  parameters  for  HCl  and  HF  discussed 
in  Section  III  are  combined  with  the  AFCRL  line  parameter  collection  evalu- 
ated in  Section  II  and  used  with  the  program  described  in  Ref,  1 to  compute 
the  radiance  from  a path  through  two  samples  of  hot  gases  characteristic  of 
possible  plume  conditions.  The  radiance  at  the  plume  boundary  is  evaluated, 
and  the  radiance  available  after  passing  through  an  absorbing  atmospheric  path 
is  computed.  The  ratio  of  spectrally  degraded  apparent  radiance  to  plume 
radiance  demonstrates  the  importance  of  including  the  effects  of  the  source 
spectrum  when  considering  the  propagation  of  radiation  through  the  atmosphere. 

The  conditions  for  the  two  gas  samples,  designated  here  plumes  A and 
B,  are  summarized  graphically  in  Figs.  21  and  22.  The  conditions  represented 
by  the  solid  curves  were  supplied  by  R.  H.  Lee"^  and  are  based  on  complicated 
plume  flow  field  modeling  studies.  These  curves  describe  the  radial  profile 
of  temperature  and  species  concentration  of  an  assumed  cylindrically  symmetric 
gas  sample  at  an  axial  location  in  the  sample  near  the  position  of  peak  infrared 
emission.  The  samples  differ  mainly  in  composition.  The  radiating  species  in 
plume  A are  CO,  CO^,  and  H^O,  whereas  those  in  plume  B and  HCl  and  HF. 
Both  plumes  are  assumed  to  be  at  a constant  pressure  of  55.3  mb. 

The  computation  of  radiance  for  inhomogeneous  samples  such  as  those 
discussed  here  is  approximated  by  a sequence  of  homogeneous  segments. 

The  values  selected  for  this  simulation  are  shown  in  Figs.  21  and  22  by  the 
horizontal  lines.  The  segments  are  not  the  same  geometric  length,  but, 
rather,  the  total  number  of  molecules  in  each  segment  is  approximately  the 


R.  H.  Lee,  private  communication  (1975). 
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different  ordinate 


Figure  22,  Temperature  and  Species  Concentrations  for  Plume  B.  The  solid  curves  denote 
the  radial  profiles  of  the  concentrations  of  HF,  HCl,  and  temperature.  The 
horizontal  line  segments  indicate  the  homogeneous  values  of  these  parameters 
used  to  simulate  the  radial  profiles.  Note  that  each  of  the  three  parameters 
has  a different  ordinate  scale. 


same.  Earlier  studies  have  indicated  that  this  provides  a better  simulation 
than  equal  geometric  length  segments.  The  numerical  values  used  in  the 
computations  are  given  in  Tables  8 and  9 for  half  of  the  assumed  symmetric 
plume . 

The  atmospheric  model  used  to  attenuate  the  radiances  consists  of  a 

path  from  20  km  to  space  at  a 60 -deg  zenith  angle  through  a mid -latitude 

3 1 

summer  model  atmosphere.  Continuum  absorption  in  the  atmosphere  was 
ignored. 

All  computations  of  spectra  were  carried  out  with  the  INHOM  program^ 
at  a resolution  of  0.01  cm  ^ for  the  entire  spectral  region  from  2 to  5 pm 
(2000-5000  cm~^).  Samples  of  the  plume  radiance  for  plumes  A and  B are 
shown  in  Figs.  23  and  24.  Apparent  radiance  spectra  were  obtained  by  mul- 
tiplying the  high- resolution  plume  radiance  spectrum  by  the  high -re solution 
atmospheric  transmittance  spectrum.  The  result  was  a spectrum  of  the 
plume  radiance  that  reached  the  top  of  the  atmosphere.  The  atmospheric 
radiance  was  not  added.  The  resulting  high-resolution  radiance  spectra 

Table  8.  Plume  A Simulation  Conditions 


Column  Densities 

Radial  Limits, 

Temperature, 

10^^  mol/cm^ 

-1.89  to  +1.89 

1403.45 

371.531 

78.866 

10.587 

1.89  to  5.023 

1165.40 

205. -i23 

47.679 

2.546 

j 5.023  to  6.601 

586.48 

64. 962 

15.553 

0.263 

t 

1 6.601  to  7.352 

279.38 

12.827 

3.074 

0.063 

7.  352  to  7.966 

227.82 

2.271 

0.627 

0 

7.966  to  8.534 

211.37 

0.341 

0.324 

U 

Table  9.  Plume  B Simulation  Conditions 


Column  Densities, 
10^^  mol/cm^ 


Radial  Limits, 
m 

Temperature, 

K 

HCl 

HF 

-1.614  to  +1.614 

1245.01 

598.402 

3642.22 

1.614  to  4.007 

923.31 

324.999 

1994.34 

4.007  to  5.017 

389.54 

80. 383 

494.963 

5.017  to  5.627 

235.37 

9.  170 

56.405 

5.627  to  6.  188 

216.70 

0.  230 

1.409 

6.  188  to  6.75 

216.70 

0 

0.001 

were  convolved  with  triangular  instrument  functions  of  two  different  widths, 

2 and  20  cm'^  FWHM,  and  are  shown  in  Figs.  25  through  34.  The  low- 
resolution  atmospheric  transmittance  spectra  are  shown  in  Figs.  25  and  26; 
the  low -re solution  plume  A radiance  spectra  are  shown  in  Figs.  27  and  28; 
the  low -resolution  apparent  plume  A radiance  spectra  are  shown  in  Figs.  29 
and  30;  the  low -resolution  plume  B radiance  spectra  are  shown  in  Figs.  31 
and  32;  the  low -resolution  apparent  plume  B radiance  spectra  are  shown  in 
Figs.  33  and  34.  Figures  31  through  34  have  the  same  vertical  scale  in 
order  to  emphasize  the  reduction  in  average  radiance  of  plume  B compared 
with  that  of  plume  A because  of  the  lower  density  of  lines.  Thirty-two  times 
the  values  of  the  20-cm”^  resolution  radiance  is  also  plotted  for  assistance 
in  determining  quantitative  values  if  desired. 

The  dependence  of  the  effective  atmospheric  transmittance  spectrum 
on  the  source  spectrum  is  emphasized  in  Figs.  35  and  36.  The  solid  curve 
represents  tlie  20-cm  resolution  average  transmittance  identical  to  the 
curves  in  Figs.  25  and  26.  The  dashed  and  broken  curves  were  obtained  by 
calculating  the  effective  average  transmittance 
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Figure  36.  Dependence  of  Apparent  Transmittance  on  Source  Spectrum 
Between  3200  and  4000  cm"^.  The  solid  curve  denotes  the 
average  transmittance  at  20  cm~^  resolution  from  Fig.  26. 
The  dashed  curve  is  the  effective  average  transmiitance  for 
Plume  A,  which  is  the  ratio  of  the  20  cm"*  resolution  curve 
of  Fig.  30  to  that  of  Fig.  28.  The  broken  curve  denotes  the 
effective  average  transmittance  for  Plums  B,  which  is  the 
ratio  of  the  20  cm"*  resolution  spectrum  from  Fig.  34  to 
that  of  Fig.  32. 


L(v)  T(v)  dv 


L(v)  dv 


where  L(v)  is  the  source  radiance,  and  t(v)  is  the  atmospheric  transmittance 
both  at  infinite  resolution.  is  essentially  the  ratio  of  the  20 -cm  resolu- 
tion appa  ent  radiance  curve  to  the  appropriate  20 -cm  ^ resolution  plume 
radiance. 

The  accuracy  of  spectra  presented  in  this  section  is  limited  by  several 
approximations  in  addition  to  the  validity  of  the  plume  models  for  any  real 
system.  As  demonstrated  in  Section  II,  the  AFCRL  line  atlas  is  not  com- 
pletely adequate  for  modeling  hot  gases.  DF  and  DCl  have  been  omitted 
from  the  plume  B computations.  The  pressures  for  the  gas  samples  are 
on  the  border  line  of  the  validity  of  the  pres  sure -boradened  Lorentz  line 
shape.  A Voigt  line  shape  was  not  used  because  of  the  additional  computer 
expense. 
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APPENDIX.  LINE  PARAMETERS  FOR  HCr  , HCr  , AND  HF 


The  line  parameters  calculated  as  described  in  the  text  for  T = 296  K 
are  listed  in  Table  A-1  in  order  of  increasing  frequency.  Listed  are  all 
lines  of  the  0-1,  1-2,  2-3,  3-4,  and  0-2  bands  of  HCl^^,  HCl^^,  and  HF  with 
strengths  such  that  at  3000  K,  a 1-m  path  of  pure  gas  at  atmospheric  pres- 
sure would  nave  an  absorption  of  1 percent  or  greater.  The  format  of  the 
computer  listing  of  each  line  is  defined  in  Table  A-2. 


T.»ble  A-1.  Line  Parameters  for  and  HF 


\m:lH 

\V4:ni 

1 310. '■•14 
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134ft,  35«; 
1346.714 
1377,134 
1 371,131 

1 3*R,4ft7 

1338. «57 
1314.43  3 
1316.730 
1471.672 
1471. ^77 
1411.774 
14*  1.617 
14?2.677 
147?, ^17 
1471,143 
1471,386 
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Table  A-2.  Format  of  Spectroscopic  Line  Atlas  Listing 


Characters 

FORTRAN 

Format 

Parameter 

Units 

1 - 10 

FiO.3 

Line  center  frequency  - 

-1 

cm 

11  - 20 

ElOiS 

Line  strength  at  296  K 

cm  / 2 

(molecule  • cm^  ) 

21  - 25 

F5.3 

Line  -width  parameter 
at  273  K 

cm  ~ ^ / atmosphere 

26  - 55 

F10.3 

Energy  level  of  lower 
state  in  transition 

-1 

cm 

36  - 44 

Blank 

45 

11 

Vibrational  quantum 
number  of  upper  state 

None 

46  - 56 

Blank 

57 

11 

Vibrational  quantum 
number  of  lower  state 

None 

58  - 65 

Blanlc 

66 

A1 

Rotational  branch 

None 

67  - 69 

13 

Rotational  quantum 
number  of  lower  state 

None 

70,71 

Blank 

72,73 

12 

”74'*  for  date 

None 

74,75 

Blank 

76,77 

12 

Isotope 

35 

15  HCl 
17  HCl^^ 
19  HF 

None 

78,79 

Blank 

80 

11 

Molecule 
3 HCl 

. 7HF- 

None 
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'Im*  l.almrilnry  0|i«ran«int  «•(  'I ha  Aarr>a|ta«a  CorpitritUin  It  ronA'iclinit 
ta|iarliiM>ntal  and  Ihaiirallcal  invaatl|a(l»na  naraattry  lor  lha  avaluttlon  and 
afipll)  aiiiin  Ilf  adanilflc  advanrat  lu  naw  mllllary  conrapit  and  tytlamt,  Var^ 
aalUlly  and  ftaalhlllly  hava  liaan  davalopad  lo  a hl||h  dagraa  hy  lha  tab>iralnry 
paraiifitiKl  In  daallng  with  lha  many  prohlarnt  anfoiiniarad  In  lha  ntllun't  rapidly 
davaliiptn  i apart  and  mlaalla  ayalama,  Kaparllaa  In  lha  lalaal  aelanllfic  daval> 
upmanla  la  vllal  In  lha  accnmpllahmani  nf  laaka  ralalad  In  lhaaa  prohlama,  Tha 
lahnral'irlaa  lhal  rnnlrlhula  In  Ihia  raaaareh  arai 

Aarnphyaira  liatinrainf vi  haunch  and  raaniry  aarndynamica,  haal  lrana» 
far,  riTanlry  phyalea,  ehamlral  VIntllct,  alruclural  machanica,  fll|hl  dynamica, 
Bimnapharlio  iKillnilon,  and  hl||h«powar  gaa  laaara, 

f hamialrv  and  Phvaira  Lahnratoryi  Almoapharic  raaclluna  and  almoa* 
pharic  '/iiilc a,  chaml cal  raaclliina  in  potlutad  aimnapharaa,  ehamlcal  raacllona 
nf  aacitad  apaclaa  In  rorlial  pluriiaa,  chamlaal  lharmndynamic a,  plaama  and 
|aaar*iMdiicad  raaetiona,  taaar  chamlalry,  proputalon  rhamlalry,  apaaa  vacuum 
and  radlaiinn  affaeta  on  malarlala,  luhrlaallon  and  aiirfara  ohanomana,  pholo' 
aantillva  malarlala  and  aananra,  high  praciainn  taaar  ranging,  and  lha  appll' 
ralinn  of  phyalea  and  chamlalry  to  prohlama  of  law  anforcamanl  and  blomadlclna, 

r-l? <:.} W Mh Klaeiromagnallc  lhaory,  davleaa,  and 
prnpagailon  phanumtna.  Including  plaama  ataclromagnallcai  quantum  alacironlea, 
laaara,  and  alaclro«oallcai  eommunlaallon  aclanrat,  applied  alacironlea,  taml> 
cnndu'  Hng,  tuparconoucllng,  and  cryalal  davica  phyalea,  optical  and  acoutileal 
Imaging!  aimoapharh'  polluiinni  mllllmalap  wave  and  far>lnfrarad  iaehnotogy, 

Malarlala  gclancag  Lihtiralnfyi  DavaWqrmanI  of  naw  rnalarlalat  malal 
rnalrla  ( nrnpnailaa  andnaw  forrria  of  carhoni  laal  and  avtiuallon  of  graphite 
and  laramlca  In  raantryi  apacaerafi  malarlala  and  tiacironic  componania  In 
nuclaar  waapona  anvironmanli  appllrallon  of  fraelura  machanica  lo  alraaa  ror> 
rnainn  and  fallgua>lnduead  fraeluraa  In  alruclural  malalt, 

finaca  Phyalea  hahoraloryi  Almoapharic  and  lonoapharle  phyalea,  radla- 
lion  fr,)m  ilia  atmoaphara,  danally  and  eompnalilon  of  lha  atrnoaphara,  aurorae 
and  alrglowi  magnalnapharlc  phyalea,  coamie  raya,  ganarallon  and  oropagallon 
nf  plaama  wavaa  In  lha  magnaloapharai  aniar  phyalea,  aludlaa  of  tolar  magnallc 
flaldti  apace  aaironomy,  n>ray  aalrnnomyi  lha  affacia  of  nuelaar  anploalona, 
magnallc  alurma,  and  aolar  artivliy  on  lha  aarlh'a  almoaphara,  lonoaphara,  and 
magnaloapharai  lha  affacli  of  optical,  eltelromagnallc,  and  parllculala  radia* 
Ilona  In  apaca  on  apaca  lyalarna. 
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